Replica Field Theory of the Dynamical Transition in Glassy Systems 
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The critical behaviour of the dynamical transition of glassy system is controlled by a Replica Sym- 
metric action with n = 1 replicas. The most divergent diagrams in the loop expansion correspond 
at all orders to the solutions of a stochastic equation leading to perturbative dimensional reduction. 
The theory describe accurately numerical simulations of mean-field models. 



More that twenty years ago it was recognized pj] that 
a certain class of mean-field spin-glass models (techically 
speaking, those where replica symmetry is broken at one 
step) displays the same dynamical behaviour predicted 
by the Mode-Coupling Theory (MCT) ofglasses [1]. This 
has motivated a great deal of research |3| and many be- 
lieve that there is an intrinsic analogy between these two 
classes of systems. 

Lowering the temperature, these models display two 
transitions. At the dynamical transition temperature Td 
the paramagnetic equilibrium state abruptly splits into 
a number of states that is exponentially large in the size 
of the system. Correspondingly the equilibrium dynam- 
ics displays the well known MCT phenomenology. At a 
static transition temperature < Td the number of equi- 
librium states is no longer exponential and this is strongly 
reminiscent of the entropy crisis which is supposed to 
happen at the Kauzmann temperature in glasses. Be- 
tween Td and T^, the entropy crisis supplemented by nu- 
cleation arguments is used to explain the observed super- 
Arrhenius slowing down of the dynamics and there are 
currently many efforts to verify this scenario. 

On the other hand MCT offers a good description of 
the early stages of the dynamical slowing down and in 
this Letter we discuss finite-size and finite-dimensional 
effects at Td, i.e. at the Mode-Coupling temperature. 
The dynamical transition is characterized by a diverg- 
ing correlation length and therefore it is natural 
to ask which are the mean- field critical exponents, what 
is the upper critical dimension above which they are 
valid, and how they are renormalized below. It is worth 
noticing that we should rather speak of pseudo-critical 
exponents because, due to activated effects, the dynami- 
cal transition is bound to disappear or rather to become 



a cross-over in finite dimension, much as in more conven- 
tional metastable phenomena. 

In spin-glass models the local variables are usually 
Ising spins Si = ±1 and the relevant order parameter is 
the overlap between different equilibrium configurations 
q(cr, r) = The same order parameter can 

be studied for liquids in a lattice gas representation. It 
is well known that in mean-field models the dynamical 
transition can be located through a static potential V{q) 
defined as the average free-energy cost of imposing that 
the system stays at an overlap q from a fixed equilibrium 
configuration The potential has a minimum at the 
equilibrium value of the overlap (q = in the absence of 
external fields), but below Td a secondary minimum ap- 
pears at a higher value, q = qd, in the mean-field theory 
of the problem. Using the replica method the potential 
V{q) can be expressed as a field theory of an action of 
a replicated order parameter that is an n x n overlap 
matrix qat with qaa — (to be continued analytically 
to n — ^ 0). The paramagnetic solution is qab = down 
to Ts, however at Td a solution appears with qab — qd 
inside an m x m block with m = 1 (analytic continua- 
tion) d, 0| ■ We want to study the loop expansion around 
the latter solution which encodes the presence of an ex- 
ponential number of states. Such an expansion can be 
simplified noticing that only the modes with diverging 
Gaussian fluctuations are relevant for critical behaviour, 
while all the others can be integrated out. Since Td is 
essentially a spinodal point for the m = 1 solution, it is 
natural to assume that the only critical variables qab are 
those inside the m x m block. We thus reach the con- 
clusion that the relevant field-theory of the problem in 
the Ginzhurg- Landau sense is a cubic Replica-Symmetric 
(RS) field theory with n — 1 replicas: 
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where (pab is the difference between the actual value of 
the order parameter qab and its mean-field value q^^^ 
and 0aa = 0. In the above expression we have retained 



only the cubic terms relevant for critical behaviour. In 
the critical region we have mi oc {T^^ — T)^/^, where 
T^^^ is the dynamical temperature in the mean-field ap- 
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proximation. The corresponding dressed propagators can 
be associated to the fohowing physical quantities: 



Gi{x,y) = {4>ab{x)(l)ab{y))c = (Sj;Sy)2 - (s^)2 (s^)2 (2) 



G2{x,y) = {(f)ab{x)(j)ac{y)) c = (Sx-Sy) {s^) (Sy) - (Sx)^ (Sy)^ 

(3) 



Gz{x,y) = {<pab{x)(l)cd{y))c = (sx)^(sy)2 - {sxy {SyY 

(4) 

In the r.h.s. of the above expressions the overhne means 
averages over the quenched disorder and over the ex- 
ponential number of states in which the paramagnetic 
state splits at Td, while the various angle brackets mean 
thermal averages computed inside the same equilibrium 
state. We can also consider connected correlations w.r.t. 
the thermal noise inside the same state, by making linear 
combinations of Gi , G2 and Gs that cancel the last term 
in eqs.dlHll), e.g. 



GsG{x,y) 
Gth{x,y) 



Gl - 2G2 + G3 = {{SxSy) ~ {Sx){Sy))^ , 

G1-G2, G,ix,y) = Gi-G3 (5) 



While these connected correlations describe fluctuations 
inside a given state, Gi, G2 and G3 yield fluctuations 
between different states, and will be called disconnected 
correlations in the following. 

We note that the above theory can be used also to de- 
scribe systems without quenched disorder, notably struc- 
tural glasses. The key requirement is that there is an ex- 
ponential number of equilibrium states and in this case 
the overline have to be interpreted just as an average over 
them. This so-called self-induced disorder is supposed to 
be the bridge between spin-glasses and glasses and it is 
the reason why the replica method can be successfully 
applied to glasses [8|, iS] . 

At the Gaussian level, fluctuations are controlled by 
the three eigenvalues of the quadratic part of C: replicon 
Tij, longitudinal and anomalous rA 10]. Degenera- 
cies occurs at special values of n: n = — >■ = r^, 
n — I r^ — rn, n — 2 ^ ta — r^. We have found 
that, as a consequence of the degeneracy at n ~ 1 be- 
tween the replicon and longitudinal eigenvalues, a double 
pole appears in the bare propagators. Switching to mo- 
mentum representation the leading divergent behaviour 
is given by: 
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where rn — r^ — 2mi and a = lim„_).i(r/j — ri)/(n— 1) — 
—2m2 — 2m3. Note that r^ is not critical and remains 
non-zero. The presence of an unexpected double pole is 
similar to what happens in the Random Field Ising Model 



(RFIM) [Tl[. Another feature of the above expressions 
that resembles the RFIM is the fact that any connected 
correlation diverges instead with a single pole. It is well 
known that the perturbative loop expansion of the RFIM 
is the same of a stochastic equation [12l] : quite surpris- 
ingly we have found that this property is also shared by 
the RS field theory with n = 1. We stress that in the lit- 
erature similar loop expansions are usually limited to the 
first few orders [ifl [HI for general values of n because for 
each diagram one has to perform a complex summation 
over replica indices. Therefore it is remarkable that in 
the case n = 1 one controls the loop expansion at all or- 
ders. The most divergent diagrams in the loop expansion 
of the theory ([T]) corresponds at all orders to the solution 
of the following cubic equation in presence of a quenched 
Gaussian random field h{x) [19] 



- A(t) + rii(j)-^ij(j)^ = h , 
with [h{x)h{y)]h — 2a6{x 



(6) 
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There is an unique cubic constant w = 0^2 — wi and 
means average over the random field. The precise mean- 
ing of the equivalence is that the most divergent diagrams 
in the loop expansion of Gi (x, y) , G2 (x, y) , G3 (a;, y) coin- 
cide to all orders with the loop expansion of [(j){x)(j){y)]h 
where (p is solution to eq. (O. The correspondence also 
holds at the level of the less divergent connected correla- 
tions functions and we have: 
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[^h]h , Gth = 2[^h]h , Gg = 3[^h]h (7) 



The stochastic equation ^ leads to a simple physical in- 
terpretation: critical behavior is controlled by the random 
field fluctuations from state-to- state. The diagrammatic 
analisys shows that the upper critical dimension is c?„ = 8 
and not the naive expectation 6. Furthermore the map- 
ping to the stochastic equation implies a perturbative 
dimensional reduction |12l |. suggesting that the critical 
exponents are the same of the pure model in dimension 
D — 2. In the RFIM dimensional reduction does not hold 
because of non perturbative effects [l3| while it does hold 
for branched polymers [l^ . There is no general recipe to 
know if it holds or not and we will not further comment 
on this point. 

At the mean-field level the theory predicts critical ex- 
ponents different from those of a standard cubic theory. 
Remarkably the predictions of the latter have been found 
to disagree with numerical simulations in a recent study 
[Tsj . In mean-field it is natural to study the critical be- 
havior of integrated quantities diverging at the critical 
point, like the susceptibilities x, = yG,{x,y), 
where N is the system size. As usual the loop expan- 
sion can be recast formally in order to deal with diver- 
gences at criticality, and the result is that disconnected 
susceptibilities diverge as N^/"^, while connected suscep- 
tibilities diverge as N^^^. However the prefactors are 
expressed as series with all positive coefficients, and are 
not resummable. This can be seen also noticing that the 
stochastic equation ([5]) does not admit a real solution 



3 



for a sufficiently negative field meaning that the averages 
[. ■ are not well defined beyond perturbation theory. 
This is precisely what we were expecting because we are 
dealing with metastable states that are intrinsically ill- 
defined at the critical temperature. On the other hand, 
dynamics is always well defined, and allow us to access 
the critical region (e.g. in numerical simulations) to test 
the above critical exponents. 

We have studied numerically the dynamics at Td, using 
a mean- field model of N Ising spins Si = ±1, interact- 
ing by 3-spin couplings {Jijk = ±1), randomly chosen 
such that each spin participates exactly to z = 8 inter- 
actions. Starting from an equilibrated configuration, we 
measured the overlap C{t) = X]il=i Si{0)si{t)/N between 
the initial configuration and the configuration at time 
t. Below Td the initial condition determines the state in 
which the system will remain trapped along the dynam- 
ics. The fluctuations of C{t) are called X4{t) in the glass 
literature [3]. However, according to our results it is cru- 
cial to distinguish between disconnected and connected 
fluctuations; so we write Xii't) — Xhet(t) + Xth(t) with 

XHedt) ^ N (Jew -WW)') , (8) 

XtH{t) ^ N {cm - {Cit)y . (9) 

In the t — !■ oo limit, we have that Xhet{t) X2 = 
N-'E.,yG2{x,y), xthit) ^ xth = N-^j:^yGthix,y) 

and X4{t) Xi = ^~^J2x y^'^i-'^^v)- '^^i^ connection 
allow us to test the above scaling predictions through 
dynamical measurements. 

Numerically we work on the Nishimori line, such that 
the starting configuration with all Si(0) = 1 is an equi- 
librium configuration [l6l |. Therefore for each sample 
we consider just one of the exponentially many typical 
states and averages between different states are obtained 
changing the sample. Averages are computed by a num- 
ber of samples (states) Ns such that NNs — 3 • 10^ 
and with 2 real replicas per sample evolving with dif- 
ferent thermal noises. The dynamics exhibits the MCT 
phenomenology. In particular above Td the average cor- 
relation Cav{t) = {C{t)) relaxes to zero in a two-step 
fashion. On the time scale tp of the (3 regime Cav{t) re- 
mains around a plateau value Cp, while one the larger 
time scale ta of the a regime it decays to zero. Both the 
two time-scales diverge at Td and the system remains in 
the state selected by the initial configuration. The an- 
alytical solution of the model ^ gives Td = 1.3420(5) 
and Cp = 0.750(5). Even at T = Td, due to finite size ef- 
fects, Cav{t) eventually decays to zero although on time 
scales diverging with the system size. In present con- 
text we do not want to discuss the behaviour of these 
time scales and we find convenient to consider reparam- 
eterized quantities XhetiCav) and Xth{Cav) [Hi- We may 
distinguish three regions. 

The perturbative region corresponds to Cav > Cp. It is 
reached in times of order 0(1), such that fluctuations are 
themself 0(1) and admit a regular expansion in powers of 
1/A^. Using rather natural finite-time-scaling arguments 
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FIG. 1: (color online) Rescaled disconnected susceptibility 
N~^^'^Xhet(Cav) (top) and rescaled connected susceptibility 
N-''^%h{Cav) (bottom). The curves for different system 
sizes cross at the analytically known plateau value Cp — 0.75. 
The insets show the same data rescaled horizontally according 
to the scaling variable x = N~^^'^(Cav — Cp), leading to a very 
good data collapse. 



one can match these dynamical perturbative series with 
the perturbative series of the statics generated by action 
(HI) . It can be argued that the coefficients of the various 
terms in the expansions diverge with the same power of 
SC = {Cav — Cp) of the corresponding static quantity. 
In particular the most divergent terms at any order are 
given by 

^ oo 

X.edC, + SC) = ^E(^^, (10) 
1 rl 

X..(Cp + 5C) ^ ^E(^. (11) 

The scaling region corresponds to values ofCav{t) near 
Cp. This region is explored on a time scale diverging with 
the system size and we expect diverging fluctuations. Ac- 
cording to the above perturbative series the scaling re- 
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FIG. 2: Normalized disconnected susceptibility X4 in the a 
regime, Cav < Cp. 



gion corresponds to SC = 0{N ^/^), with the following 
scaling laws 

XhetiCp + N-^/^x) = N^'^fueAx), (12) 
Xth{Cp + N-"^x) = N^'^ftH{x), (13) 

where the two scaling functions fhet{x) and fth{x) go to 
zero for x — +oo and diverge for x — )■ — oo. The above 
perturbative series provide their asymptotic expansion 
for X — ^ cjo: at leading order fhet(x) ~ a;~^ and fth{x) ^ 
x^^. The numerical data are in very good agreement 
with the expected behaviour, see Fig. [TJ 

The a regime corresponds to < Cav < Cp and 



we can not use any perturbative information here. In 
this region our numerical data are close to the form 
X4(Cai)) = NCav{Cp — Cav), scc Fig. [21 which would hold 
if, for any initial configuration and thermal noise, the re- 
laxation consisted in a sharp jump from the plateau value 
Cp to the value for uncorrelated configurations, C = 0. 

It is worth noticing that it is far from trivial that 
the exponents predicted in perturbation theory can be 
actually observed in the critical region. Indeed if non- 
perturbative effects were also present on the same time 
scale they would wash out the perturbative results. Our 
results suggest that non-perturbative effects appear in- 
stead on a larger time scale. The time scale of the crit- 
ical region diverges with the system size as iV*^ and 
a reasonable possibility is that the time scale of the a 
regime scales like A^i^ + ie ^ where a and b are the dynam- 
ical MCT exponents of the /3 regime. The results of the 
static theory can be used to safely infer other properties 
of critical dynamics in the early /? regime. In particular, 
ioi T = T ~ Td > 0, we have a /3-like regime for short 
time, t < tfj = 0(t~^/^"), where the following scaling 
laws should hold 



Xhet{t) ^ - ghetit/tfj) , xth{t) ^ —=gth{t/tii) . (14) 
r 

Correspondingly, we have X4(0 Xhet(i) <x. t^"" and 
Xth{t) oc at finite times. It could be that one is able 
by means of matching arguments and numerical observa- 
tion to access also the late /3 regime and the a regime. 
We choose however not to discuss this point because we 
think that a satisfactory understanding of the a regime 
shall include an analitycal treatment of non-perturbative 
effects. 
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